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ABSTRACT 

The 4 Ms Chandra Deep Field-South (CDF-S) and other deep X-ray surveys 
have been highly effective at selecting active galactic nuclei (AGN). However, 
cosmologically distant low-luminosity AGX (LLAGN) have remained a challenge 
to identify due to significant; contribution from the host galaxy. We identify 
long-term X-ray variability ( —month years, observed frame) in 20 of 92 CDF-S 
galaxies spanning redshifts z cs 0.0S — 1.02 that do not meet other AGN selection 
criteria. We show that the observed variability cannot be explained by X-ray 
binary populations or ultraluminous X-ray sources, so the variability is most 
likely caused by accretion onto a supermassive black hole. The variable galaxies 
are not heavily obscured in general, with a stacked effective power-law photon 
index pfT 5lK k as 1.93±0.13. and are ther efore likely LLAGN. The LLAGN tend to 
lie a factor of «6 80 below the extrapolat cd linear variability-luminosity relation 
measured for luminous AGN. This may be explained by their lower accretion 
rates. Variability-independent black-hole mass and accretion-rate estimates for 


department of Astronomy & Astrophysics. 525 Dave;,' Lab, Tlic Pennsylvania State University, Univer- 
sity Park, PA 16802, USA 

2 Dipartimento di Seien/e Fisiehe, Universita' Ffcderico II di Napoli, Via Cinthia, 80126 Napoli, Italy 
'department of Physics, University of Durham, South Road. Durham, DH1 3LE, UK 

1 Pont; ices Universidad Catolica de Chile, Departamento de Astronomfa y Astroffsica, Casilia 306, San- 
tiago 22. Chile 

-Space Science Institute, 4750 Walnut Street. Suite 205, Boulder, Colorado 80301 
'The Johns Hopkins University, Homewood Campus, Baltimore, MD 21218, USA 
7 NA$A Goddard Space Flight Centre, Code 662, Greenbelt, MD 20771, USA 
8 Harvard-Smithsonian Center for Astrophysics, 60 Garden St. Cambridge, MA 02138 USA 
9 Department of Physics, University of North Texas, Denton, TX 76203, USA 
10 llniversita di Bologna, Via Ramani 1, Bologna, Italy 


variable galaxies show that they sample a significantly different black-hole mass 
accretion rate space, with masses a factor of 2.4 lower and accretion rates a 
factor of 22.5 lower than variable luminous AGN at the same mlshift. We find 
that an empirical model based on a universal broken power-law PSD function, 
where the break frequency depends on SMBH mass and accretion rate, roughly 
reproduces the shape, but not the normalization, of the variability-luminosity 
trends measured for variable galaxies and more luminous AGN. 

Subject headings: galaxies: active X-rays: galaxies 


1. Introduction 

Observations show that all nearby galaxies with a massive bulge component host super- 
massive black holes (SMBHs) (c.g., Ferrarese & Ford 2005; Giiltekin et al. 2009). SMBHs 
accreting near the Eddington limit (L/L Kdd ~ 0.1 .1) are visible as luminous active galactic 
nuclei (AGN) that often outshine their host galaxies. Models of AGN lifetime, constrained 
by observed Eddington ratio distributions, suggest that SMBH growth is dominated by this 
luminous phase, lasting — a few x 10 s years (c.g.. Marconi et al. 2004: Shankar et: al. 2004; 
Hopkins & Hernquist 2009). 

Observations including constraints on the sizes of ionized “bubbles" around quasars 
(c.g., Jakolisen et al. 2003; Gonealves et al. 2008) and the lengths of relativistic j»t.s and radio 
lobes (e.g., Scheuer 1995; Blundell et al. 1999) suggest that the episodic lifetime of luminous 
activity is similar to the total lifetime, implying that a SMBH is triggered to the luminous 
AGN phase no more than a few times. SMBHs therefore spend significant amounts of time in 
quiescent or low-activity phases, which may contribute up to —20% of overall SMBH growth 
(Hopkins k Hernquist 2009). A SMBH accreting at lower rates (L/Lmu 0.1) will appear 
as a low-luminosity AGN (LLAGN). LLAGN share several properties with more luminous 
AGN. including similar X-ray spectral shapes (e.g., Younes et; al. 2011) and similar radio- 
loud fractions and luminosity-dependent., intrinsic X-ray to optical flux ratios (e.g., Maoz 
2007), A more complete census of LLAGN is important for understanding SMBH accret ion 
history, but the relative significance of the host galaxy in LLAGN makes a full census of 
accretion activity a challenge. 

Deep X-ray surveys have been effective at selecting a wide variety of AGN, including lu- 
minous, unobscurcd AGN as well as faint and/or obscured AGN (c.g., see Brandt k Hasingci 
2005, for a review). X-ray selection criteria usually include cuts on X-ray luminosity and 
X-ray spectral shape. Multi-wavelength data further aid X-ray selection by allowing seine* 
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lion via, for example, excess X-ray emission compared to what is expected from optical flux 
(e.g., Hornschemeier et al. 2003) or radio luminosity (e.g., Alexander ct al. 2005). 

The above methods have been successful in selecting a wide variety of AGN, but nev- 
ertheless miss certain populations, such as very heavily obscured AGN and LLAGN (e.g., 
Bauer ct al. 2004; Alexander et al. 2005; Lehmer et al. 2008). While heavily obscured AGN 
can often be selected in the IR (e.g., Houck et al. 2005; Martinez-Sansigre et al. 2005: Alexan- 
der et al. 2008), the spectral energy distributions (SEDs) of LLAGN are likely dominated by 
the host galaxy in other bands. Even in X-rays. X-ray binaries, ultraluminous X-ray sources, 
and hot gas will provide significant contributions to the overall power output. Simulated 
Chandra observations of nearby low-luminosity Seyfert nuclei artificially shifted to 2 ~ 0.3 
show that LLAGN would exhibit X-ray luminosities, spectral shapes, and X-rav-to-optical 
flux ratios consistent with those of normal or optically bright/X-ray faint galaxies (Peterson 
et al. 2006). By relying on such criteria, deep X-ray surveys may be underestimating AGN 
fractions. 

Variability potentially provides a useful indicator of whether an extragalactic X-ray 
source, classified as a galaxy by other means, harbors an AGN. Variability is a defining 
characteristic of AGN and has long been used as an AGN selection technique (e.g., van den 
Bergh et al. 1973). Numerous studies have used optical variability to select AGN from deep 
surveys such as the 1 Ms CDF-S, the Subaru/XMM-Newton Deep Field, and the GOODS 
North and South Fields (e.g., Trcvese et al. 2008; Morokuma 2008; Villforth et al. 2010; 
Sarajedini et al. 2011, respectively). Spectroscopic observations of the 1 Ms CDF-S (Boutsia 
et al. 2009) found that 17 of 27 optical variability-selected objects were broad-line AGN; 9 
(5) AGN would have been missed if selected by color (X-ray selection). 

Similarly, UV variability has been used successfully to identify LLAGN in galaxies with 
low-ionization nuclear emission-line regions (LINERs). LINERs have been found in the 
nuclei of a large fraction of nearby galaxies (e.g.. Ho et al. 1997; Kauffman et id. 2003), but 
these regions could be ionized by cither massive star clusters or low accretion-rate AGN. 
HST imaging has found that ~25% of LINERs are associated with compact (< few pc) UV 
sources (Maoz et al. 1995; Barth et al. 1998). A study of LINERs with compact nuclear UV 
sources found significant variability in 15 of 17, indicating the presence of an AGN (Maoz et 
al. 2005). 

Deep X-ray surveys are able to detect variability in moderate-luminosity /high-redshift 
AGN (e.g., Almaini et al. 2000; Paoliilo et al, 2004; Papadakis et al. 2008b). The 4 Ms 
Chandra Deep Field-South (Xue et al. 2011), the deepest X-ray survey to date, allows a 
preliminary classification of AGN on the basis of several observed quantities (see §2 for 
details). This paper utilizes X-ray variability techniques to search for AGN missed by these 


criteria. With 4 Ms of exposure time spanning 10.8 years for 466 good-quality sources (see 

§3), variability can be detected in sources with time-averaged fluxes as faint as F 0 , 5 8tev » 

5 x 10~ 17 ergs cm -* s' 1 . We use a cosmology with Ho = 70.4 km s' 1 Mpc ’, Q fl , = 0.272, 

and S! a = 0.728 (e.g., Komatsu et al. 2011). 


2. Overview of the 4 Ms CDF-S Catalog 

The details of the 4 Ms CDF-S source catalog are available in Xue et al. (2011); we 
provide a brief summary here. The 4 Ms CDF-S, constructed from 54 Chandra observations 

over 10.8 years, covers an area of 464.5 aremin 2 and reaches highest sensitivities of Fa . 5 2tev ~ 

9.1 x 10-“ ergs cm" 2 s' 1 and F 2 _ 8tev » 5.5 x 10 17 ergs cm ' 2 s' 1 , with multi-wavelength 
coverage in more than 40 bands from the radio to the UV. Source candidates are detected 
using a 10 -5 false-positive probability threshold in WAVDETECT (Freeman et al. 2002) and 
are then pruned using a binomial no-source probability (see Appendix A2 of Weisskopf et al. 
2(1)7) P < 0.004 to obtain a more conservative list of 740 main-catalog sources, all of which 
are consistent with being point sources. Source extraction and photometry were conducted 
with ACIS EXTRACT (AE; Broos et al. 2010). AE models Chandra' s High Resolution Mirror 
Assembly point spread function (PSF) using the MARX ray-tracing simulator. 11 The PSF 
model is used to generate a polygonal extraction region for each sounc that approximates 
the «90% encircled energy fraction (EF,F) contour of a local PSF measured at 1.497 keV. For 
<6% of the candidates, the sources are crowded (i.e. , the polygonal source regions overlap) 

and smaller extraction regions that are as large as possible without overlapping (40 75% 

EEF) are used. The background is calculated from regions that subtract the contribution 
from the source of interest and its neighboring sources; the regions are typically a factor 
«16 larger than the source-extraction region. AE merges the individual observations to 
estimate aperture-corrected, background-subtracted counts and the ler (asymmmetric) upper 
and lower statistical errors (Gehrels 1986). In this paper, we will use the standard X-ray 
photometric bands: 0.5 2 keV (soft). 2 8 keV (hard), and 0.5 -8 keV (full). 

Though most sources have a relatively small number of counts (median net counts » 
77). a rough estimate of source spectral shape can be made by relating the band ratio (i.e., 
the ratio of the count rates in the 2 8 keV and 0.5-2 keV bands) to an effective power- 
law photon index, F c | f (F„ oc v~ a = i/" ,,+1 ). For low-count sources where r,, ff cannot be 
determined reliably, r,, (f is set to 1.4, the stacked (co-added) sped rum of all sources in (he 
CDF-S (Tozzi ct al. 2001: Xue et al. 2011), which is consistent with the unresolved spectrum 


"MARX is available at http://space.vnit.edu/CXC/MARX/index.html 
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of the cosmic X-ray background (Hickox k Markevitdi 2006). 

Of 740 X-ray sources, 716 (96.8%) contain matches in at least one of seven optical/near- 
infrared/radio (ONIR) catalogs: (1) the ESO 2.2 m WFI A-band catalog (Giavalisco et 
al. 2004), (2) the GOODS-S Hubble Space Telescope (HST) version r2.0z z-band catalog 
(Giavalisco et al. 2004), (3) the GEMS HST z-band catalog (Caldwell et al. 2008), (4) the 
GOODS-S MUSIC catalog (Grazian et al. 2006). (5) the MUSYC K- band catalog (Taylor 
et al. 2009), (6) the SIMPLE Spitzer/IRAC 3.6pm catalog (Damon et al. 2011), and (7) the 
VLA 1.4 GHz radio catalog (Miller et al. 2008). 

Of 716 X-ray sources with multi-wavelength identifications, 419 (58.5%) have spectro- 
scopic redshift measurements, collected from Le Ffcvre et al. (2004), Szokolv et al. (2004), 
Zheng et al. (2004), Mignoli et al. (2005), Ravikumar et al. (2007), Vanzella et al. (2008), 
Popesso et al. (2009), Tfeister et al. (2009), Balestra et al. (2010), and Silverman et al. 
(2010). A total of 343 (81.9%) of the 419 spectroscopic redshift measurements are “secure", 
in that they are measured at > 95% confidence levels with multiple secure spectral features. 
668 (93.3%) sources have high-quality, accurate (|Az|/(l + z) ss 6.5%) photometric-redshift 
measurements from at least one of three photometric-redshift catalogs: Cardamone et al. 
(2010), Luo et al. (2010), and Rafferty et al. (2011). The positions of primary ONIR coun- 
terparts were Cross-matched with the photometric-redshift catalogs using a matching radius 
of 0.5", resulting in a fake-match probability of <1%. Subsequent spectroscopic observations 
published in the Arizona CDF-S Environment Survey (ACES: Cooper et al. 2012) catalog 
show with a blind test that errors on the photometric redshifts are <1%. 

The 4 Ms CDF-S X-ray sources were classified as AGN by the following criteria: 

• High luminosities: Lo.o-skev > 3 x 10 42 ergs s -1 , where the rest-frame luminosity has 
been corrected for Galactic and intrinsic absorption. 

• Hard spectra: A source with rW < 1.0 is identified as a heavily obscured AGN. 

• High X-ray-to-optical flux ratios: log (Fx/Fr) > - 1, where Fx = Fo.s-skev, Fo. 5 _ 2 kev, 
or Fi-gkev and /■% is the /7-band flux. 

• Excess X-ray emission compared to that expected from star' formation: Lo.s-stev > 
3x(8.9x 10 17 Ln) (Alexander et al. 2005), where Lh is the 1.4 GHz monochromatic luminosity 
in W Hz" 1 . 

• An indication of broad emission lines in the optical spectrum. 

Stars were classified by cross-matching X-ray sources (using the ONIR counterpart po- 
sitions) with (1) the spectroscopically identified stars in Szokolv et al. (2004). Mignoli et al. 
(2005), and Silverman et al. (2010); (2) the likely stars with stellarity indices greater than 
0.7 in the GEMS HST catalog (Caldwell et al. 2008); and (3) the likely stars with best-fit 
stellar templates in the MUSYC photometric-redshift catalog (Cardamone et al. 2010), using 


a matching radius of 0.5". 

X-ray sources not identified as an AGN or a star were classified as galaxies. 

Rest-frame 0.5 8 keV luminosities are calculated for all sources. For AGN, which make 
up the vast majority of CDF-S sources, the luminosity is corrected for Galactic and intrinsic 
absorption. AGN not detected in the full band have an upper limit on the X-ray luminosity 
based on the 3cr Poisson error on the counts. For galaxies, the intrinsic-absorption correction 
and A'-corrections may not be appropriate. Of 92 CDF-S galaxies with “good quality” ob- 
servations (sec §3), 78 tire not detected in the hard band and have poorly determined photon 
indices. All galaxies are detected in the soft-band. We perform a stacking analysis for galax- 
ies with < 150 net counts (87 of 92) following the procedure described in §5.1. The resulting 
average = 1.90±0.08 gives ias-skcv/Aos 2 kev - 2.35. We calculate the rest-frame 
0.5-8 keV luminosity for CDF-S galaxies as Lo.s-skeV 7 4ir (P L x 2.35 x /o.s- »*v(l + z) r “**- 2 . 
Since galaxy X-ray emission is typically unabsorbed, we do not apply any correction for 
intrinsic absorption. 


3. Testing for X-ray Variability 

We perform two quality cuts before conducting variability analysis. First, we exclude 
the catalog sources with off-axis angles greater than 8' to ensure that sources will have 
sufficient coverage (>50 of 54 observations) t hroughout the various CDF-S pointings. To 
ensure accurate variability measurements, wo also require that each source has at least 20 
net counts in the 0.5 8.0 keV band (i.e., at least 5 counts on average per epoch, as defined 
below). These quality cuts result in a total sample of 466 CDF-S sources: 369 classified as 
AGN, 92 classified as galaxies, and 5 classified as stars in the Xue et al. (201 1) catalog. All 
of these AGN and galaxies have measured spectroscopic or photometric redshifts. The 92 
sources classified as galaxies may nevertheless contain an AGN. as would be indicated by 
X-ray variability, and make up the sample investigated in this paper. 

Wc divide the CDF-S observations into four epochs, each containing ~1 Ms of integra- 
tion: 2000 (943.1 ks), 2007 (DC7.7 ks), 2010a (March May; 1015.5 ks), and 2010b (May July; 
944.9 ks). As in Xue et al. (2011), we merged observations within each epoch and, for a 
given source position from the CDF-S catalog, measured the source and background counts 
and flux over three observed-frame energy bands: 0.5 8 keV, 0.5 2 keV, and 2 8 keV. A 
source is considered variable if the variability observed between observations is greater than 
that expected from Poisson statistics, with a probability threshold of 5% that the variability 
is due to noise alone. (The choice of probability threshold is discussed further below.) To 
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check whether a source can be considered variable, we calculate the quantity: 


A’ 




fa - !‘f 


( 1 ) 


where N = 4 is the number of epochs, ./■, is the photon flux (background-subtracted counts 
with units of cm 2 s _1 ) in a given epoch, /( is the mean photon flux over all epochs, and aj 
is the error squared on the photon flux for the i‘ h epoch. The photon flux is calculated by 
dividing the full-band (0.5-8.0 keV) net counts by the exposure time and the mean effective 
area across the source aperture. The Gehrels (1986) approximation gives the error on the 
net counts, which is propagated to obtain the error on the photon flux. Since this error is 
significantly asymmetric for low-count (< 15 counts) sources, we average the upward and 
downward error bars for these objects to obtain the average error cr,. (The same method is 
applied in the Monte Carlo simulations below.) 


For large photon fluxes, the X 2 statistic follows a X 2 distribution, and any source with 
X 2 > 7.82 (for 3 degrees of freedom) has a probability P X 2 < 0.05 (i.e., 95% confidence level) 
that the variability is due to random noise. However, at low count rates, the error on the 
photon flux is not Gaussian. Since errors in the low-count regime are larger than expected 
from a Gaussian distribution, the resulting A' 2 statistic is smaller than expected and docs 
not follow the X 2 distribution (see Fig. 2 of Pnolillo et al 2004). 


We therefore constructed a Monte Carlo simulation to determine the distribution the 
X 2 statistic should follow for each source, similar to the procedure followed by Paolillo et 
al. (2004). We first scaled the total observed source and background counts for each source, 
obtained from the full 4 Ms observation, to the exposure time and effective area for a given 
epoch. This procedure generates the source and background counts expected in each epoch 
if the source and background were constant over time, and it accounts for fluctuations in the 
background that will affect low-count sources. To simulate the variance expected from noise, 
Poisson distributions were defined using the expected source and background counts as the 
mean values. We then simulated 1,000 observations of each source by repeatedly drawing 
the expected counts from the Poisson distributions. For each simulated observation, we cal- 
culated the photon flux for four epochs and calculated A' 2 as defined above. Asymmetric 
errors on the source and background counts are obtained from Gehrels (1986) and are prop- 
agated to get the error on the photon flux. The observed X 2 can then be compared to the 
simulated distribution to determine the probability P x t that the observed variability is due 
to Poisson noise. Spurious sources of variability are negligible, since effective exposure maps 
are calculated separately for each observation, taking into account issues such as vignetting, 
CCD gaps, bad pixels, bad columns, and Chandra s spatial- and time-dependent quantum 
efficiency degradation. 


While using a more conservative P crit — 1% on our dataset would result in fewer false 
positives, it would also eliminate a similar number of truly varying sources. For example, in 
the sample of 92 galaxies, P„ A ~ 5% results in 20 variable sources (sec §4), of which 4.6 are 
expected to be false positives. Reduc ing P C!il to 1% results in 13 variable sources, of which 
0.9 is expected to be a false positive. So while the more conservative critical value eliminates 
«4 false positives, it also eliminates » 3 4 truly varying sources. The P x t values listed in 

Table 2 can be used to screen the sources further as desired. 


4. Galaxies with AGN-like Variability 

Of the 369 CDF-S sources classified as AGN that, meet both the off-axis angle (<9 < 8') 
and count (total net counts > 20) requirements, 50.1% exhibit significant flux variability 
(P x j < 0.05) on —month year timescales. For the 178 AGN with more than 100 counts, 
74.2% exhibit significant variability. The basic diagnostic plot in Figure 1 shows the rise in 
the AGN variable fraction with total net counts. The plot demonstrates that, given sufficient 
counts to detect it, variability is a near-ubiquitous property of faint, X-rav selected AGN. 
even in the case of significant obscuration: a 70% of the CDF-S sample consists of obscured 
AGN (Xue et al. 2011), and 47.5'% (51.5%) of obscured (unobscured ) AGN arc significantly 
variable. The variable fractions are consistent with the results from t.he 1 Ms CDF-S (Paolillo 
et al. 2004), although obtained with different temporal sampling and down to much fainter 
fluxes. 

In sources classified as galaxies, the variable fraction is significant at low counts and 
equals that of AGN at higher counts (Figure 1). Table 1 describes the columns of Table 2, 
which lists the attributes of the 20 variable and 72 non-variable galaxies. (Variability prop- 
erties of CDF-S AGN will be covered in the forthcoming Paolillo et al., in preparation.) The 
luminosity- redshift distribution of variable and non-variable galaxies is shown in Figure 2; 
almost all lie below z ~ 1. Spectroscopic redshifts are available for 18 of 20 variable galaxies 
and for 61 of 72 non-variable galaxies. Photometric redshifts are available for the remaining 
sources. Six example light curves (background-subtracted count, rates in the observed-frame 
0.5-8 keV band vs. M.JD) representative of the sample as a whole are shown in Figure 3. 

The X-ray luminosity distributions of all CDF-S galaxies and those exhibiting significant, 
variability are shown in Figure 4. A K-S test shows that the two samples are consistent with 
being drawn from the same parent population (P Hl! ~ 57%), and the variable fraction does 

not show a significant dependence on X-ray luminosity below 7r. -, ~ 10 4:! ergs a '. The 

possibility that AGN-relatcd variability may go undetected in galaxies is discussed in §4.1. 
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We briefly compare AGN selection based on variability to the following selection meth- 
ods: (1) X-ray luminosity cuts, (2) the X-ray-to-optical flux ratio, and (3) excess X-ray 
emission compared to that expected from star formation, based on the radio luminosity. 
Figure 5 shows the fraction of variable sources vs. X-ray luminosity. Below Lo-s-skev = 
10 42 ergs s' 1 , a luminosity cut often used for AGN selection in X-ray surveys, the variable 
fraction remains significant, at 20 - 30%. Of 04 CDF-S galaxies with i-as-swv < 10 42 ergs s _l , 
17 (ss27%) are variable. 

AGN selection via the X-ray-to-optical flux ratio is demonstrated by the /?-band mag- 
nitude vs. X-ray flux plane in Figure 0 (cf. Figure 16 in Xue et al. 2011). Sources classified 
as AGN. galaxies, and stars in Xue et al. (2011) are marked as small red circles, larger black 
circles and blue stars, respectively. Variable sources are marked with filled symbols. The op- 
tically bright., X-ray faint region, typically thought to exclude AGN (OBXF; Fx/Fr < —2; 
Hornschemeier et al. 2003) contains 27 galaxies. Of these, 6 (22%) are variable, with X-ray 
luminosities spanning log 7o y sov -- 39.7-41.4. 

AGN may also be selected based on a comparison between X-ray and radio luminosities 
(Xue et al. 2011). The radio luminosity can be used to predict the X-ray luminosity in star- 
forming galaxies (e.g., Alexander et al. 2005, and references therein), so sources with excess 
X-ray emissioif may tx* classified as AGN. Of the 17 CDF-S galaxies with radio detections 
(none of which have excess X-ray emission), three (18%) are variable. These objects may 
have excess radio emission due to strong radio cores rather than star formation. 


4.1. Measuring Variability Strength 

Due to the generally limited photon statistics of the CDF-S galaxy sample, most vari- 
able sources must be strongly variable to be detected. Significantly variable galaxies show 
maximum-to-minimum flux ratios Ti max / m in « 1.5 9.3 (median = 4.1) over the observed 10.8- 
year time frame. The smallest max-to-min ratio (1.5) was measured for 033246.77-274212.7 
(XID = 616), a source with >500 counts. For most galaxies, total net cotints are too low 
(< 100) to detect variability below a factor of 2-3. 

To address whether the variable galaxy max-to-min ratios sample the average AGN 
population or only the highly variable “tip of the iceberg.” we ran a Monte Carlo simula- 
tion. The procedure assumed that the entire galaxy population is significantly variable and 
simulated the variability expected over 10.8 years of observation. Following the procedure 
in Vaughan et al. (2003), we used the Timmer k Koenig (1995) algorithm to simulate 5,000 
light curves based on the mean and variance of the flux for each of the 92 CDF-S galax- 


ies. For non-variable galaxies, the measured variance represents noise in the measured light 
curve, which gives an upper limit to the variability that could be present. The algorithm 
produces a random, continuously sampled light curve from a given power spectral density 
(PSD) function, which we assumed to be described by a broken power-law, where the break 
frequency depends on SMBH mass and accretion rate (McHardy et al. 2006). Since the 
break frequency lies outside the range of timescales sampled for most variable galaxies (sec 
§6.3), wc simplify the model to a power-law (P(f) a f 0 ) with index ft = 1, as is typical 

for the low-frequency (long-timescale) PSDs for nearby Seyferts (e.g., Vaughan et al. 2003). 
The time units of the light curves are determined by the minimum and maximum timescales 
input into the simulation; the light curve durations are adjusted according to each source’s 
redshift. 

Wc resampled/rebinned the light, curve using the CDF-S observing pattern and add 
Poisson noise to the simulated light curve to account for measurement error. The full sim- 
ulated light curves was made five times longer than the sampled region in order to produce 
variation on timescales much longer than those sampled by the data. This reproduces the 
effect that very long- timescale (low-frequency) variations have on variability measured over 
shorter timescales (i.e., “red noise leak"). 

We compare the simulated distribution of median max-to-min ratios, where the median 
is calculated over 5,000 trials for each galaxy, to the observed distribution for variable galaxies 
in Figure 7. A K-S test shows that these two populations have a 0.2% chance of being drawn 
from the same parent population. Note that the simulated distribution (black histogram) 
illustrates an upper limit to the variability that could be present in the CDF-S, resulting in 
a lower limit on the detectable fraction of sources. The CDF-S detects at least »18% of 
sources with max-to-min ratios > 2 and at least «57% of sources with max-to-min ratios > 
4. A significant fraction of non-variable galaxies may still host an AGN, but the variability 
may remain undetected due to low counts. 

For the galaxies exhibiting significant variability, wc calculated the normalized excess 
variance (e.g., Vaughan et al. 2003), which measures how strongly a source varies in excess 
of the measurement error. The excess variance is the integral of a source’s PSD function 
over a given frequency range, which is defined by the light curve’s duration (10.8 years, 
observed-frame) and minimum bin size (4.0 Ms, 12 observed-frame). 


12 While each CDF-S epoch totals —1 Ms in integration time, the timescale sampled is significantly longer 
due to the spread in Chandru observations, of which the shortest is 4.0 Ms for the 2007 epoch {Luo et a!. 
2008). 
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where a„ ri is the average of the asymmetric upward and downward measurement errors. 
Using the upward (downward) error would overestimate (underestimate) the errors. Zero 
excess variance (<rjj x , = 0) would indicate that the observed count fluctuation is entirely 
consistent with noise rather than due to intrinsic source variability; due to statistical fluc- 
tuations. the excess variance may also be negative in this case. The ir 2 xs values are listed in 
Table 2. Note that the variability amplitude is calculated for observed-frame energy bands. 
The variable galaxies cover a redshift range z ~ 0 1, so the excess variance will be mea- 

sured over 0.5 - 8 keV at z = 0 to 1 - 16 keV at 2 = 1. Variability strength is known to 
depend on energy in some nearby Seyfert galaxies (e.g., Ark 120, MCG 6-30-15, and I Z w 1; 
Vaughan et al. 2004, Vaughan & Fabian 2004, and Gallo et al. 2007, respectively), while in 
others, variability remains nearly constant with energy (e.g.. Ton S180 and NAB 0205+024; 
Vaughan et al. 2002 and Gallo et al. 2004, respectively). In the former cases, the change in 
variability strength is small, with a < 10% difference from 0.5 to 10 keV, so the bandpass 
effects at redshifts 2 =0-1 should remain small. 


We calculate the statistical error (i.e., measurement error) on the excess variance fol- 
lowing Equation 11 of Vaughan et al. (2003). 


err(er 2 xb .) = 



( 3 ) 


The large errors on cr 2 xs (Table 2) are due to the small numbers of counts observed for most 
sources (e.g., 12 of 20 variable galaxies have < 100 net counts); four variable galaxies with 
net counts < 50 have excess variance measurements completely dominated by statistical 
uncertainty [o 2 xs < crr(o 2 xs )). Nevertheless, most variable galaxies have excess variance 
measured at the >1 a level. 

The excess variance contains additional sources of uncertainty aside from statistical 
error: (1) random scatter intrinsic to the stochastic nature of AGN variability (Vaughan et 
al. 2003) and (2) uncertainty and systematic bias due to sparse sampling of the light curve. 
The sparse pattern of CDF-S observations will lead to large uncertainty in the mean flux 
measurement, and since the measured mean will be closer to the sampled data points rather 
than the true mean, the variance measurements will tend to be underestimated (Allcvato et 
al. in preparation). We again employ a Monte Carlo simulation to model the uncertainty 
and systematic effects produced by random scatter and sparse sampling. 


We follow the same procedure described above to produce 5,000 light curves for every 
variable source. The mean and variance arc calculated after sampling each simulated light 
curve with the CDF-S observing pattern. The sampling bias can be corrected by rescaling 
the observed variance by a factor equal to the ratio between the “true” input variance (i.e., 
the observed variance used as input in the simulations) and the median output variance (i.e., 
the biased variance affected by sparse sampling): = CT 2 eM /median((T s 2 im ). The median 

output variance is calculated over all 5,000 light curves. The amount of systematic bias 
depends on the frequency and regularity of the sampling. If the sampling is regular, the 
scaling factor will approach unity as the number of samples increases; however, the scaling 
factor will remain above unity even at high sampling frequency if the sampling is irregular 
(Allcvato et al. in preparation). The slope of the PSD will also affect the sampling bias 
a steeper PSD slope (i.e., 0 = -2 instead of 0 - -1) will result in a larger bias for a 
given sampling pattern. Since the intrinsic PSD slopes are not known, we apply uniform 
corrections assuming 0 = —1: however, source- to-souree variations in PSD slope may be a 
significant source of scatter in variability measurements. 

We find that the scaling factors range from f mle « 0.03 to 4.8, with a mean 
fa 1.54 and a scatter on /., ,v. of Or rs 0.87. In —76'/ of the sources, the scaling factor is 
greater than unity, indicating that the variance is underestimated due to the Sampling bias. 
Note that the Monte Carlo PSD is normalized bv each source’s variance, which is calculated 
using the source’s light curve. Since the heavy binning could smear out high frequency 
variations, this method could result in an artificially smaller scaling factor. However, Allcvato 
et al., in preparation, find a similarly small bias (/«, i e < 2) when sampling higher frequencies 
for a wide range of S/N ratios, gap lengths and sampling patterns. More importantly, the 
scatter on the /so*, factor calculated for each source is large (« 40% 190%), so individual 

measurements, even when corrected for bias, are likely to be poor estimates of the intrinsic 
variance. Therefore, variability properties of CDF-S galaxies are best considered in ensemble 
rather than on an individual basis. Bias-corrected excess variances (<t 2 ) are listed in 

Table 2 and are used for all further analysis unless otherwise noted. 


4.2. Comparisons with XRB and ULX Variability 

The three most likely sources of X-ray variability are X-ray binary (XRB) populations, 
ultrahiminous X-ray sources (ULXs), and accreting SMRHs. tit this section, we show that 
the first two possibilities are not likely to dominate the measured galaxy variability. 

To examine the potential contribution of X-ray binary populations to variability, we 
must first determine the relative contributions of low-mass X-ray binaries (LMXB) and high- 
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mass X-ray binaries (HMXB) to the hard (2 10 keV), galaxy-wide X-ray luminosity (Lxrb). 13 
Galaxy stellar mass (M t ) scales the contribution of (older) LMXBs, and star formation rate 
(SFR) scales the contribution of (younger) HMXBs (e.g., Equation 3 of Lehmer et al. 2010). 
M t and SFR are calculated for each galaxy in Xue et al. (2010) using the optical colors and 
the UV and IR luminosities. Both the Lehmer et al. (2010) relations and the Xue et al. 
(2010) calculations adopt the same initial stellar mass function (Kroupa 2001). However, 
since Lehmer et al. (2010) and Xue et al. (2010) use different formalisms for computing 
stellar masses (Bell et al. 2003 and Zibetti et al. 2009, respectively, which differ primarily 
due to their models of star formation history), we apply a correction factor of 2.6 to the 
stellar masses from Xue et al. (2010). The SFR and AL* values for each variable galaxy are 
listed in Table 2. 


By comparing the expected LMXB and HMXB contributions to the total luminosity, we 
can determine which population ought to dominate the variability. In the variable galaxies, 
SFR ranges from 0.04 to 55 Mq yr -1 with a median of 2.2 Mq yr 1 : AL* ranges from 2 x 10 7 
to 6 x Hi 11 M. with a median of 1.2 x 1O 1O M 0 . We find that all but six of 20 variable galaxies 
are expected to have a larger HMXB contribution. Lhmxb(SFR)/Z, LMX b(A/*) ranges from 
0.01 to 75.8 with a median of 2.6 (Table 2). HMXBs are generally more variable than LMXBs 
(Gilfanov et ah 2004), so unless LMXBs dominate the X-ray output of a galaxy, we neglect 
their contribution. 


To determine the variability expected from the HMXBs, we follow the relations in §4.2.3 
of Gilfanov et al. (2004), where the variability of the HMXB population is roughly determined 
by the galaxy’s SFR. The following relations were obtained from Monte Carlo simulations 
with a power-law HMXB luminosity function with slope a = —1.6 and a cut-off luminosity 
at L cut = 2 x I0 40 ergs s~‘: 


„ 0.35l“ io for SFR < 5 AL 0 vr" 1 (4) 

Ormn.O 

~ 0.30 + g J? for SFR = 5 10 AL 0 yr -1 (5) 

rms,0 

1- for SFR > 10 M e yC 1 (6) 

Orm.fi SFR'" 


Here, o rtm ,o is the fractional rms (i.e., the square root of the excess variance) expected from 
an individual X-ray binary, which can be as large as 20-30% on ~year timescales (e.g., 


13 We measure variability in the 0.5 8 keV band, but by limiting the comparison to the 2 10 keV luminosi- 
ties, we limit the contribution of other sources of galaxy-wide X-ray emission (hot gas, supernovae, supernova 
remnants, and O-stars). which fade sharply above 2 keV and can be considered negligible. 


Gilfanov 2010), and o rms , it , is the total variability. We take o rmSi 0 = 0.3 and calculate the 
upper limit on a> rms , M . In the most extreme case, we find that the upper limit, on variability 
expected from an HMXB population is erJ KB < 0.014. 

In six galaxies, LMXBs me expected to dominate the X-ray luminosity. Four of the six 
have log AL. > 10.5 AL 0 , and therefore have largo enough stellar mass to follow a « 

Miteufr I aw (§4.4.2 of Gilfanov .et al. 2004). For the remaining two galaxies, we follow the 
trend in Figure 12 of Gilfanov et, al. (2004). The upper limit on variability expected from 
an LMXB population is < 0.02. 

We find that XRB populations cannot explain the full extent of the X-ray variability 
for variable galaxies. Figure 8 plots the distribution of for GDI' S galaxies, showing 

that all the variable galaxies exhibit variability in excess of that expected from XRBs. The 
variable galaxies have a median o™ iC<>rr /(Tx fl B ~ 42 (without the bias correction, a 

14), indicating that the contribution of an XRB population to the measured variability is 
small. 

We also consider whether one or more ultraluminous X-ray sources (ULXs) may domi- 
nate a galaxy’s X-ray output. The nature of ULXs is still debated, but most likely involves 
accretion onto massive stellar black holes (30 100 A/,.,); a few cases may involve aecretion 
onto intermediate mass black holes (100 300 Af@) or beamed emission from HI 20 M . black 
holes. Typical luminosities span Lo.s-skeV ~ 10“ 10 41 ergs s _l (e.g.. Lehmer et al. 2006: 
Swartz et al. 2011). In a survey of 1,441 X-rav point sources in 32 nearby galaxies in the 
Chandra archive, Colbert, et al. (2004) found ULXs in 19 galaxies; the contribution of one 
or more ULXs to a galaxy’s total X-ray point source luminosity ranged from 7% 87%. with 
a median contribution of 43%. 

ULXs could potentially explain nine variable galaxies with Los sk, ,v < 10 41 ergs s 1 . 
However, since ULXs tend to be associated with star-forming regions (e.g., Swartz et al. 
2004) and occur more frequently in late-type/irregular galaxies than in earlv-typc galaxies 
(e.g., Walton et al. 2011), the five variable galaxies with both Los-ski-v < 10 41 ergs s 1 and 
late-tvpe morphology (see §5.2) are more likely to host ULXs. 

One method of finding ULXs is to search for off-nuclear sources. We plot postage- 
stamp images (8" x 8") of the GOODS-S/GEMS HST V 606-band for the variable galaxies 
in Figure 9. The circle overplotted on each image has a radius 1.5 times the Chandra 
positional erro, which is calculated at the 90% significance level. As in Lehmer et al. (2006), 
X-ray sources offset from the galaxy nucleus by more than 1.5 times the positional error are 
considered off-nuclear. 


We find one marginally off-nuclear source: 033219.27 275406.7 (X1D 269). The 
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primary optical source appears to be an early-type galaxy with blending toward the galaxy 
to the lower right, suggesting a merger. Both galaxies have similar redshifts (2 = 0.960 and 
0.956, respectively); the first redshift is spectroscopic and secure (see §2) and the second is 
photometric (Xuc et al. 2010). The high X-ray luminosity (To. 5 - 8 keV « 2.2xl0' 12 ergs s' 1 ) 
suggests that a ULX is likely not the dominant source of X-ray emission from this galaxy. 

One variable galaxy, 033230.00-274405.0 (XID = 418), was previously identified as 
being off-nuclear in the 1 Ms CDF-S (Lehmer et al. 2006). With the additional data from 
the 4 Ms CDF-S, the X-ray source position has been refined (with reduced uncertainty) to 
be consistent with the galaxy’s nucleus. The off-nuclear source discussed in the previous 
paragraph (XID = 269) was not detected in the 1 Ms CDF-S, so it was not considered by 
Lehmer et al. (2006). 

We cannot rule out that a ULX may lie too close to a galaxy's nucleus to be detected as 
an off-nuclear source. However, the possibility of ULXs in most variable galaxies is mitigated 
by high X-ray luminosities and/or early-type morphology, so since we find no plausible off- 
nuelear ULXs, we conclude that ULXs are unlikely to dominate the emission from variable 
galaxies. Accretion onto a SMBH remains the best explanation of variable galaxies. 


5. Supporting Evidence for LLAGN 

We investigate the X-ray spectral shapes, the morphologies, and the optical spectral 
classifications of the variable galaxies for two purposes: (1) characterizing the variable 
galaxy population, and (2) determining whether their properties sire consistent with those 
of LLAGN. 


5.1. X-ray Spectral Shape 

As discussed in §2, Xue et al. (2011) calculate the effective photon index for each source 
based on the ratio of count rates in the 2 8 keV and 0.5 2 keV bands (Table 2). For 10 
variable galaxies detected in the soft band but not the hard band, lower limits are listed. For 
4 low-count variable galaxies detected cither in the full-band, soft-band, or both, no reliable 
effective photon index can be determined, so F„ir is set to 1.40. 

Since 14 sources have poorly determined F e if, we perform a stacking analysis on all 
variable galaxies with < 150 net counts to determine an average photon index for the sample. 
The three highest count sources (with 199.2, 275.6, and 608.9 net counts) are excluded from 
the stacking analysis since they could dominate the stacked signal, but the results remain 


the same within errors if these sources are included. 

Following the procedure in Luo et al. (201.1), we calculate the soft and hard-band counts 
in a 3" diameter aperture for each source. The background is calculated by averaging the 
counts in 1,000 randomly placed, source- free apertures within a 1 -radius circle around the 
source position. The individual source counts are summed and background is subtracted. 
Aperture corrections are applied, averaged over al] the observations weighted by exposure 
time (see Luo ct al. 2011 for details), before calculating the band ratio. 

The stacked effective photon index for the 17 relevant variable galaxies is F. . - « 
1 .93±0. 13, which is consistent with the typical photon index for local Seyfert, galaxies (F 
~ 1.8; e.g., Dadina et al. 2008) at the ler level. Including all 20 variable galaxies, F, 

~ 1.82±0.08. While absorption may still be present in some individual sources, the soft 
spectrum implied by the stacked effective photon index indicates that the average variable 
galaxy is not heavily obscured. The X-ray luminosities are therefore intrinsically low and 
indicate that variable galaxies are most likely LLAGN. 


5.2. Galaxy Morphologies 

Postage-stamp images (8" x 8") of the GOODS-S/C! EMS IIST 17606-band for CDF-S 
variable and non-variable galaxies are shown in Figures 9 and 10, respectively, with Chan- 
dra error circles overlaid. Galaxies are classified by eye as late-tvpe, early-type, irregular, 
or undetermined. Mergers are also visually classified based on blending and/or tidal tails 
between two or more galaxies. 

Since visual classifications are subjective and are particularly unreliable for distant, 
poorly resolved galaxies, we also apply the color-magnitude relation given in Bel! et al. 
(2004), where galaxies are considered part of the “red sequence” (i.e. , early-type morphology) 
if they are redder than (Mu My) = -0.31s 0.08Afv 0.51. Galaxies with bluer colors 

are considered part of the "blue cloud” (i.e., late-type morphology). Although the color- 
magnitude diagram leads to a more objective classification scheme, there are nevertheless 
significant uncertainties, both in the rest-frame magnitudes and colors, and in the definition 
of the Bell et al. (2004) relation. Moreover, the color classification is complicated by the fact 
that many of the galaxies lie in between the red sequence and blue cloud, in the so-called 
“green valley.” In cases where the color and the visually classified morphology disagree, 
we choose a final classification, preferring the visually classified morphologies in nearby, 
well-resolved galaxies, and the color classification in distant and/or poorly resolved galaxies. 
Considering only galaxies classified as early or late-type, the visual and color classifications 
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agree «54% of the time. Table 2 lists each galaxy’s morphological type as determined by 
visual classification, the rest-frame magnitude and color (from Xue et al. 2011), the color 
classification according to Bell et al. (2004), and the final classification. 

Based on the final classification from Table 2, we find that variability does not prefer one 
morphology type over the other. Variable galaxies have 40.011^% early-tvpe morphology 
and 50.01]®;®% late-type morphology, compared to 23.6 tl'l% and 51.41? |%, respectively, in 
non- variable galaxies (the errors are ler binomial errors). If we instead apply only the Bell et 
al. (2004) color classifications, then both the variable and non-variable galaxies prefer late- 
type morphologies (80.0l®’,° s % and T2.2l|;|%, respectively). The M, and SFR distributions 
for variable and non-variable galaxies show no significant difference ( Pks > 10%). 

We find that the fraction of mergers among variable and non-variable galaxies is con- 
sistent within errors (9.1tJ°j 6 % and 20.0i|;|%, respectively). 


5.3. Optical Spectroscopic Classification 

Most (18/20) variable galaxies have optical spectral observations (Szokoly et al. 2004: 
Zheng et al. 2004; Mignoli et al. 2005: Ravikumar et al. 2007; Popesso et al. 2009: Silverman 
et al. 2010), from which spectroscopic redshifts were determined. In all cases, the optical 
spectra are classified as galaxies, showing only narrow emission lines or absorption lines. 

Szokoly et al. (2004) classify objects in more detail. Of the eight variable galaxies listed 
in the Szokoly et al. (2004) catalog, two have typical galaxy spectra showing only absorption 
lines. The remaining six are classified as having low-excitation emission lines consistent with 
H II rogion-tvpc spectra. These objects would be classified as normal galaxies based on the 
optical data alone as the presence of the AGN cannot be established. However, one of these, 
033222.78 -275224.2 (XID = 312), has sufficient signal-to-noisc in the optical spectrum to 
measure line ratios. This object is classified as a LINER by Szokoly et al. (2004) via the line 
ratio diagnostics given by Ho k Sargent (1993). 


6. Galaxy vs. AGN Variability 

6.1. The Variability-Luminosity Anti-Correlation 

In this section, we investigate how galaxy variability compares to AGN variability. We 

first confirm a significant anti-correlation between excess variance and X-ray luminosity 

among AGN. as seen in previous work (e.g., Barr & Mushotzky 1986: Lawrence k Papadakis 


1993; Nandra et al. 1997; Hawkins 2000; Paolillo et al. 2004). The anti correlation is plotted 
in Figure 11 as cr 2 * X!liCorr vs. Lo.s~si«v. The excess variance includes the bias correction 
discussed in §4.1. The rest-frame X-ray luminosities arc calculated differently for AGN 
and galaxies, as described in §2. The Spearman rank correlation coefficient (p e = -0.44) 
shows the correlation is significant at 1\ as 10~ 9 * (5.9<r) for AGN only (rod ciieles). The 
correlation increases in significance to P„ as 8 x 10” 12 (6.4tr) if all variable sources above 
Lo.5-8k«v = 10 41 ergs s~' .are considered, including those classified as galaxies (black circles 
and stars). No significant correlation is found if the sample is limited to only the variable 
galaxies, which are discussed further in §6.2. Note that the correlation coefficient and best-fit 
line parameters (given below) remain consistent within errors if calculated using the excess 
variance uncorrected for sampling bias (o 2 „ s ). 

Since both variability and luminosity may depend on other parameters, such as black 
hole mass and/or accretion rate, we fit the data using sixlin.pro, an IDL program adapted 
from Isobe et al. (1990). A least-squares bisector fit to sources with Lo r, a i*v > lO 41 ergs s' 1 , 
weighted by the uncertainties in a 2 XS|COn , results in <T* Mitx> , r oc Lo 5 liv*’- The slope and 
intercept are consistent within errors if the weights are not included. Limiting the sample 
further to sources with L o s 8k „v > 10 w ergs s“ 1 results in o 2 xs . rorr a S These slopes 

are both consistent with the results of Nandra et al. (1997), where a weighted least-squares 
bisector fit results in cr 2 x , <x L/!\o(,v 22 , which is notable since the Nandra et al. (1997) data 
sampled shorter (hour day) timescales compared to the months years timescales sampled 
by the CDF-S. 

The slope presented in this paper is significantly flatter than that found in Paolillo et 
al. (2004; o 2 IS oc which included non-varying sources. We choose not to include 

upper limits for non-varying sources in our analysis because the assumptions underlying 
survival analysis, which have been successfully applied to deal with censored data in other 
astronomical situations, do not apply here because: (1) the excess variance measurements 
of most sources lie near the detection limit, (2) a large percentage of sources do not have 
detected variability («50% of AGN and «78% of galaxies), and (3) a significant percentage 
of sources with no detected variability, especially those classified as galaxies, may truly not 
be variable (i.e., o 2 xs ~ 0). By not including censored data in this paper’s measurements, 
we likely bias the measured slopes and possibly the significance of the anti-correlation. Nev- 
ertheless, the variability-luminosity anti-correlation has been observed to follow a model 
based on SMBH mass and accretion rate (e.g., Papadakis et al. 2008b), suggesting that the 
anti-correlation is real, though probably not linear. The model is discussed further in §6.3. 

We cheek for other potential biases that may affect the luminosity- variability anti- 
correlation. The log-log plot of Figure 11 has the disadvantage of “hiding" sources with 
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negative <r^„ cort values. Sources with low flux values, and hence higher scatter in d^ !Cotr 
will therefore appear to have stronger variability, since the <Tn„^. ort values scattered to neg- 
ative values will be hidden. To check for this bias, we remove high scatter sources with 
err (a^„ corr ) > 0.1, excluding all but 38 sources, and find that the <j^ <corr — Lo s-skev cor- 
relation remains significant, at. F, «s 0.9% (2.0 a). The best-fit slope and intercept remain 
consistent at the la level. 

The flux-limited nature of the CDF-S survey presents another potential bias. Since lu- 
minosity correlates with redshift, and intrinsic variability timescales decrease with redshift, 
the decrease in variability strength as luminosity increases could, in principle, simply reflect 
the fact that shorter timescales are studied at higher redshifts, and therefore exhibit less 
variability due to the red-noise nature of AGN light curves. To check for a possible redshift 
bias, we examine the crj xscorr -i 0 . 5 - 8 kev correlation using a sub-sample within a narrow red- 
shift range (0.55 < z < 0.75). This redshift range selects 10% of the total sample, covers 
luminosities from Lo.s-skev ~ Iff 41 - 5 to 10 4 ' 1 ergs s -1 , and results in minimal differences in 
rest-frame timescales. The Spearman rank correlation coefficient for this sub-sample remains 
significant at P, ss 0.8% (2.6cr): the slope and intercept me consistent at the ler level. 

The above test also addresses another potential bias due to the redshift range of the 
sample. The a'* m values listed in Table 2 measure the variability strength in the observed 
frame, so they will sample different energy bands depending on the source’s redshift (see §4. 1 
for discussion) . This could introduce bias if the variablitv amplitude changes with increasing 
energy: from 0.5 — 8 keV at z — 0 to 1 — 16 keV at z -- 1. The slope and intercept remain 
the same within the narrow redshift range tested above, suggesting that any such bias does 
not have a significant effect on the sample. 


6.2. Suppressed Variability in LLAGN 

Variable sources with luminosities less than £o.r,-8keV = 10 41 ergs s -1 tend to fall signif- 
icantly below the extrapolated linear relation by factors of ss6 80 (median factor of «24). 
indicating a drop in variability relative to the linear relation on long timescales for LLAGN. 
This “suppressed” variability can be shown to be intrinsic to properties of AGN variability 
rather than due to dilution by unrelated XRB populations. 

In §4.2, we showed that galaxy variability cannot be attributed solely to XRB popula- 
tions. We now check whether the XRB contribution could nevertheless dilute the observed 
variability by estimating how much XRBs are expected to contribute to the total hard (2 
10 keV), galaxy-wide X-ray luminosity (Lx rb; see §4.2: Lehmer et al. 2010). 


We compare L X rb to the measured, intrinsic X-ray luminosity. First, Lo.s-stev is con- 
verted to Lj-ioibVi using the intrinsic photon index of F - 1.8 adopted by Xuc et al. (2011) 
for AGN, and the stacked photon index T slark « 1 .!) for galaxies (see §2). Figure 12 plots the 
measured 2 10 keV luminosity against that expected from an XRB population for variable 
galaxies (red stars) and non-variable galaxies (green squares). For reference, we plot in the 
same figure 369 CDF-S AGN (open orange circles), 32 local galaxies (Colbert et al. 2004), 
and 20 local luminous infrared galaxies, which arc likely to be actively star-forming (Lehmer 
et al. 2010). The red solid line shows unity, while the dotted lines mark the dispersion 
observed in Lehmer et al. (2010). 

A K-S test shows no significant difference between the relative XRB contribution in 
variable and non-variable galaxies. The CDF-S galaxies on average tend to lie above unity, 
with a median L^~ wfccv/ L xrb « 4.6. The X-ray excess in the non-variable galaxy population 
is perhaps not surprising given the flux-limited nature of the CDF-S survey, which is more 
likely to detect the objects at the bright end of the galaxy X-ray luminosity function, and 
may result in a high percentage of “contamination" by AGN. 

The median X-ray excess for variable galaxies (La-utkcv.wt/ lok. v.xnn ~ 9.2) suggests 
that XRBs contribute - 1 VA of the 2 TO keV luminosity for t lie average: variable galaxy. For 
six variable galaxies, the total X-ray emission is consistent with that expected from XRB 
emission within the scatter of the Lehmer et al. (2010) relation (see Fig. 12), suggesting that 
dilution may be more significant in these sources. Three of these sources have variability 
consistent with the linear variability-luminosity relation, while three have suppressed vari- 
ability (filled black stars in Fig. 11). Dilution by XRB variability may therefore play a role 
in suppressed variability, but cannot fully explain the extent to which o;;,, SiCon . is suppressed 
at low luminosities. 

An alternative possible explanation for the suppressed variability at low luminosities 
is a change in accretion structure. Ptak et al. (1998) found a similar drop in variability 
strength below L^-iokeV ~ 2xl0 41 ergs s' 1 in a sample of LLAGN and LINERs observed 
with ASCA on variability timescales of less than a day. The authors hypothesized that a 
radiatively inefficient, accretion flow (RIAF, e.g.. Yuan k Narayan 2604) could be responsible 
for suppressed short-timescale variability at low luminosities due to the larger extent of the 
X-ray source. This scenario would not obviously explain the reduced variability on —month 
year timescales seen here. RIAF models also predict a hard X-ray photon index due to the 
lack of an optically thick accretion disk, which provides the soft X-ray photons. The stacked 
X-ray photon index for variable galaxies (r^nn, ~ 1.93:1-0.13; §5.1) is inconsistent with this 
prediction. 

Studies since Ptak ct al. (1998) have found evidence both for “suppressed" variability 



21 


22 


in LLAGN (e.g., Ptak ct al. 2004; Markowitz k Uttley 2005; Papadakis et al. 2008a) and 
against it (c.g., Binder et al. 2009; Pian ct al. 2010). Similarly, objects such as narrow-line 
Seyfert 1 (NLS1) galaxies and the dwarf Scvfert N’GC 4395 (Mbh » 3.6 x 1O 5 M 0 ) exhibit 
“excess” vairability for their luminosity (Boiler et al. 1996; Iwasawa ct. al. 2010). However, 
when plotting variability against mass instead of luminosity (e.g., Papadakis et al. 2008a; 
Miniutti et al. 2009), such discrepancies disappear, with residual differences possibly due to 
varying accretion rates (e.g., McHardy et al. 2004; Markowitz k Uttley 2005). 

To investigate the role of SMBH mass (Mon) and accretion rate (normalized by the 
Eddington rate; rh — - M/M&u), we have obtained rough estimates for all variable sources. 
Masses and Eddington ratios for all galaxies, variable and non-variable, are listed in Table 
2. The SMBH masses are estimated via the scaling relation between Mbh and absolute 
A'-band magnitude (Graham 2007); 

log— T 1 = — 0.37(±0.04)(JVjftf + 24) + 8.29(±0.08) (7) 

Mg 

The total absolute rest-frame if -band magnitudes are derived from SED fitting in Xue et 
al. (2010) with a random scatter of <0.3 mag. An X-ray luminosity-dependent correction 
factor (Equation 1 of Vasudevan et al. 2009) corrects for nuclear emission. We assume that 
the host galaxy is bulge-dominated, a valid assumption for most AGN (e.g., Kauffman et 
al. 2003; Grogin et al. 2005). Several variable galaxies, however, are not bulge-dominated 
(nine are late-type; see Table 2), so 1 heir black hole masses may be overestimated. We 
apply a luminosity-dependent bolometric correction (Ka-iokev) to estimate the bolometric 
luminosity (Marconi et al. 2004) and calculate the Eddington ratio (Lm/Lbm — [^ 2 -iofceV 
jL 2 iokevl/[l-25 x 10 : « {Mbh/Mq)]). Note, that tire Marconi et al. (2004) correction is 
calculated explicitly for lokcv > 10 42 ergs s -1 ; we extrapolate this relation down to the 
lower luminosities of the variable galaxy sample. 

Both the Mbh and Am estimation techniques are known to have large dispersions. 
The Mhh-Mk scaling relation has a total scatter of 0.33 dex, and additional uncertainty 
will come from the luminosity-dependent correction for nuclear emission, which is based on 
template SEDs (Vasudevan et al. 2009). In addition, the assumption that all the variable 
galaxies are bulge-dominated will produce additional scatter. 14 The bolometric correction, 
too, has large scatter due to the intrinsic dispersion in the SED. The uncertainty in L w 
due to SED dispersion is ~20% for luminous AGN (Elvis et al. 1994; Richards et al. 2006). 
There is some debate regarding the similarity between the SEDs of LLAGN and luminous 


H Thf: scaling relation in Kormendv <1: Gebhardt (2001). for example, has a much larger scatter of 0.56 
dex largely because of poor bulge/disc separation (Graham 2007). 


AGN (e.g., §1; Ho 1999 and Ho 2002 vs. Maoz et al. 2007); however, the luminous AGN 
SED dispersion likely serves as a lower limit to dispersion in LLAGN SEDs. 

Variable galaxies tend to have lower accretion rates ((rh) = 4 x 10 4 ) and masses 
({Mbh) = 2.6 x 10 7 M 0 ) than variable AGN ((in) >9x 10 -3 ; (Mbh) = 6.2 x 1O 7 M 0 ), 
where we have limited the AGN sample to e < 1 for purposes of comparison. A K-S test 
shows that the differences in the m and Mbh distributions are significant: P K s « 1.7x10 5 
and P ks * 0.002, respectively. The properties of variable AGN and galaxies arc consistent 
with the range of estimates made by Babic et al. (2007) for X-ray selected, z < 1 AGN in 
the 1 Ms CDF-S, which span m ~ 10" 5 1 (median ss 0.001) and Mbh ~ Vf - 10 10 A/ . 

(median « 10 s M 0 ). 

Unlike previous studies (O’Neill et al. 2005; Papadakis et al. 2008a; Miniutti et al. 2009), 
we find no significant (anti-)correlation between or 2 ,, and M B h, most likely due to the large 
scatter in o 2 ,, and M B h measurements, combined with the. narrow range of masses probed 
(since most sources lie between M mt ~ 10 ? 10® A/ a ). It is not. surprising, therefore, that we 
find a significant anti-correlation between (T 2 r8 and in (P KS » 1.6 x 10 n , 6.3o), plotted in 
Figure 13, which is likely an artefact of the cr 2 ,,- Lo.s-.8keV anti-correlation. Nevertheless, it is 
interesting that in Figure 13 the variable galaxies (black circles) connect smoothly with the 
more luminous AGN (red circles), with no discrepancy in variability strength. This suggests 
that the factor of 22.5 difference in m may explain the “suppressed" variability in variable 
galaxies. 


6.3. Comparing the Variability-Luminosity Relation with Empirical Models 

A number of recent studies (e.g., McHardy et al. 2004; Papadakis 2004; O’Neill et al. 
2005; Papadakis et al. 2008b) have shown that X-ray variability may be determined by a 
combination of Mbh and rh, explaining the observed variability-luminosity relation. AGN 
light curves appear to be universally described by a broken power-law PSD function, where 

the break frequency depends on mass and accretion rate: - 0.029j7»h ( A/ «« / 10 6 A/..,) ’, 

where r/ is the accretion efficiency', assumed to be r? — 0.1 (McHardy et al. 2006). The 
excess variance is equivalent to the integral of the PSD between the minimum and maximum 
frequencies sampled by the data, so as long as the break frequency falls within this range, 
the excess variance at a given luminosity will decrease' with increasing M B h and increase 
with increasing m. 

We compare the Mbh and m values estimated in §6.2 to the variability measured for the 
variable AGN and galaxy populations. Assuming a universal PSD function, we derive the 
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variability-luminosity relations expected for the CDF-S sampling pattern, given a range of 
SMBH masses and accretion rates (c.g., Papadakis et al. 2008b, and references therein). The 
bolometric luminosity, which is calculated directly from Mbh and m, is converted to X-ray 
luminosity via a bolometric correction (Marconi et al. 2004). Where the break frequency lies 
outside the timescales sampled by the data (depending on the combination of m and Mbh)- 
the excess variance will remain constant. 

The model variability-luminosity relations are plotted in Figure 14 for the average ac- 
cretion rates for variable (ril — 9 x 10~ 3 ; dashed line) and non-variable (m = 4 x 10 -4 ; 
dash-dotted line) galaxies covering Man — 10 5 10*0 A/ 0 (low to high Los-skev)- Both vari- 
able populations are limited to z < 1 for comparison, and the relations are calculated at z 
= 0.5 and z — 0.7, the median redshifts for variable galaxies and AGN, respectively; us- 
ing model redshifts z — 0 or 1 resulted in negligible changes. The shape of the observed 
variability-luminosity anti-correlation, including the plateau at low luminosity, is roughly 
reproduced by the model. The plateau occurs when the break frequency lies outside the 
range of timescales sampled by the data. Unfortunately, the long timescales (especially the 
long minimum timescale of 4 Ms; see §4.1) and the large scatter prevent the models from dis- 
tinguishing between the significantly different accretion rates estimated for variable galaxies 
and AGN. 

Note that most sources are more variable than predicted by the model. We note some 
possible sources of bias: (1) The normalization of the PSD function used to calculate the 
models is based on a small sample of nearby AGN (Papadakis 2004) and may therefore not 
be representative of CDF-S AGN out to z» 1. (2) Both the models and the bias correction 
applied to the measured excess variance (§4.1) depend on a universal broken power-law PSD 
function, but the slopes of the PSD may vary between individual sources, and some sources 
may even have a second break at shorter frequencies (e.g., McHardy et al. 2007). (3) Peculiar 
sources, such as NLSls, exhibit higher variability 15 , so the variable sample may suffer from 
selection effects due to the higher sensitivity to more variable sources (§4.1). 


7. Conclusions and Future Work 

We have investigated the variability properties of sources classified as galaxies in the 
4 Ms CDF-S by dividing the observations into «1 Ms epochs. We find the following main 


I5 A notable exception is the NLS1 galaxy Aku 561 , which is the only AGN observed to date to have a 
second, low-frequency break in its power spectrum, leading to lower variability on long timescales (McHardy 
et al. 2007). 


results: 

1. Using a \ 2 -based Monte Carlo simulation appropriate for low-count sources, 20 of 
92 galaxies are selected as variable, as well as 185 of 369 AGN. Variability is effective in 
selecting AGN that do not meet other AGN selection criteria, such as cuts based on X-ray 
luminosity, X-ray-to-optical flux ratio, or excess X-ray emission compared to radio emission 
expected from star-forming galaxies (sec §3). If all 20 variable galaxies are accreting SMBH, 
then the commonly used selection criteria employed by Xue et al. (2011) miss 20 of 389 AGN 
(«5.1%), and 19 of 54 LLAGN («s35.2%; where LLAGN are arbitrarily defined as having 
io.s-skev < 10 42 ergs s ' 1 ). Due to the low counts measured for most galaxies, it is likely that 
some non-variable galaxies may still host an AGN but not exhibit detectable variability, so 
these fractions serve as lower limits. 

2. We measure variability strength with the normalized excess variance, accounting 
for measurement error, red-noise scatter, and bias due to the CDF-S sampling pattern. 
Comparing the bias-corrected excess variance to that expected from XRB populations, we 
find that XRBs cannot explain galaxy variability. 

3. The possibility of energetically significant ULXs in variable galaxies is mitigated by 
high X-ray luminosities and/or earlv-type morphology, in addition, we find no plausible 

off-nuclear ULXs. 

4. Galaxy X-ray variability is most likely associated with accretion onto a relatively 
unobscured SMBH. Though some absorption may be present, the individual and stacked 
photon indices show no indication of heavy obscuration and are consistent with the typical 
AGN spectral shape (r stltt k ~ 1.93±0.13). 

5. We confirm a significant anti-correlation between excess variance and X-ray luminos- 
ity, and find the slope and intercept to be consistent with what is reported in the literature 
for shorter timescales. We show that the sampling bias induced by the CDF-S observing 
pattern does not affect, the anti-correlation slope or significance. 

6. Low-luminosity sources (Los-stev < 10 4 ' ergs s *, largely variable galaxies) have 

“suppressed” variability compared to the extrapolated linear relation between excess variance 
and X-ray luminosity. This may be explained by their lower accretion rates. Variable galaxies 
sample a different mass-accretion rate space than the rest of the AGN population, with 
significantly lower accretion rates ((m) = 4x 10 1 ) and masses ({Mbh) 2.6 x 10 7 Af,) 
than variable AGN ((rh) = 9 x 10” 3 ; {Mbh) ~ 6.2 x lO 7 Af 0 ) in the same redshift range 
(z < 1). 

7. We find that an empirical model based on a universal broken power-law PSD function, 
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where the break frequency depends on SMBH mass and accretion rate, roughly reproduces 
the shape of the variability-luminosity anti-correlation. However, the normalization of the 
model is low compared to the data, suggesting systematic bias either within the data (e.g., 
selection effects) or the model (e.g., incorrect normalization for LLAGN at z < 1). 

We have shown that the variability measured by deep X-ray surveys is an effective 
technique for selecting cosmologicallv distant LLAGN. Extending the CDF-S to even longer 
exposures would enable detection of variability in both fainter and less variable sources, thus 
allowing better characterization of the properties and abundance of the LLAGN population. 
Within the 4 Ms CDF-S, smaller time bins could be used to search for variability and better 
characterize the variability presented here. 

In addition, follow-up optical spectroscopy is necessary to characterize the LLAGN 
population presented here. While Szokoly et al. (2004) were able to classify one object as a 
LINER, the spectra for the other sources did not have sufficient signal-to-noise to measure 
line ratios, so these objects would require either deeper optical spectroscopy or stacking of 
the currently available spectra. 
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Table 1. Overview of Columns for Properties of CDF-S Galaxies 


Column Description 


1 Sequence number in the CDF-S catalog {Xue et al. 2011) (i.e., XID) 

2 CDF-S name 

3 Redshift 

4 Method of redshift measurement (s — spectroscopic, p = photometric) 

5 Upper limit flag for net counts 

6 Total net counts in the 0.5 8.0 keV band or 3cr upper limit 

7 log Io. 5 _. 8 k«v (ergs s“ l ) 

8-10 Effective photon index (r e ff) a and corresponding errors 

11 Variability statistic (A' 2 ) 

12 Probability that X 2 statistic is due to chance (Px*) 

13 Maximum-to- minimum flux ratio (2? max / m »•„) 

14-15 Excess variance (<r 2 xs )Wd corresponding error 

16-17 Bias-corrected excess variance (cr 2 xs corr . see §4.1} and corresponding error 

18 SFR (Af® vr" 1 ) 

19 M* (10 50 M g ) 

20 Tumxb/^lmxb 

21 O'XRB 

22 Visual galaxy classification 

23 Mv 

24 Mu - Mv 

25 Color galaxy classification 

26 Final galaxy classification 

27 log Mbh (Me) 

28 log M/Mtcdd 


a The effective photon index is calculated from the band ratio. For sources detected in 
the soft band but not the hard band. F f , tf is a lower limit. When the counts are too low 
to determine reliably the photon index from the band ratio, F is set to 1.4. 

b Errors on excess variance are calculated according to Equation (11) from Vaughan et 
al. (2003). The excess variance is not significant for non- variable galaxies, but the errors 
can be used to calculate an upper limit. 


Table 2. Properties of CDF-S Galaxies 


XID 

CDF-S name 

z A 

Net 

counts 

log 

f-'O.S-BkflV 

Frit 

X‘ 

h 



120 

033206.40-274728.6 

1.02 

32.7 

41.62 

1.4 

2.8 

0.0296 

4.1 0.22 ±0.24 

0.56±0.63 

154 

033209.79 274442.7 

0.08 

46.9 

39.39 

>0.9 

2.9 

0.0304 

4.1 0. 32.3.0.21 

0.851:0.75 

162 

033210.72-274234.9 

0.42 

111.7 

41.16 

i nH’-l 
JU -0.3 

5.8 

0.0002 

• 3.4 0.22:1:0.11 

0.601.0.49 

223 

033215.80-275324.7 

0.67 

<43.7 

41.51 

>1.0 

5.4 

0.0034 

8.1 0.56:1-0. 29 

1.67:11.34 

233 

033216.76-274328.2 

0.52 

61.0 

41.26 

>1.0 

4.4 

0.0050 

3.5 0.43:10.21 

1.16:10.97 


Note. Table 2 is presented in its entirety in the electronic edition. A portion is shown here for guidance regarding 
its form and content. The full table contains 28 columns of information for the 20 variable and 72 non- variable CDF-S 
galaxies. 

a For 18 of 20 variable galaxies and 61 of 72 non-variable galaxies, the redshift is measured spectroscopically and is 
“secure” (see §2) . The remaining galaxies have photometric redshifts. 
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Fig. 1 . The fraction of variable AGN (red circles) and galaxies (black squares) vs. net 
counts, where the background-subtracted (net) counts are measured in the 0.5 8.0 keV band 
over 4 Ms. The binomial errors on the fraction at the 90% confidence level are calculated 
from Cameron (2011). The error bar on the net counts represents the bin size. The points 
are slightly offset for clarity. 


Timescale (years) 

10.8 7.2 5.4 4.3 



Fig. 2. The 0.5 8 keV luminosity vs. redshift. distribution for CDF-S galaxies meeting the 
quality criteria of §3 (open circles). Galaxies with significant variability are marked as red 

stars. The CDF-S on- axis flux limit for 20 net counts (F 0 . 5 g^y «s 4.7 x 10 17 ergs cm 2 s ' ) 

is plotted as a dotted line. The upper s-axis shows the maximum rest-frame timescale in 
years, where the maximum observed-frame timescale of the CDF 8 is 10.8 years. 
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Fig. 3. Light curves for 6 representative CDF-S sources classified as galaxies that exhibit 
significant variability. Full-band (0.5 8 keV, observed-frame) counts have been grouped into 
the four epochs and asymmetric errors on the count rates were calculated via Gehrels (1986). 
The mean count rate is overplotted as a dashed line, and the error on the mean is shown as 
the grey shaded area. Each plot lists the P* 2 values and source names. 
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Fig. 4. The distribution of X-ray luminosity (Lo. 5 - 8 k«v) for all CDF-S sources classified as 
galaxies (dark grey histogram) and for those exhibiting significant variability (red hist ogram). 
The vertical lines show the median values for all galaxies (black, dotted) and for variable 
galaxies (red, dashed). A K-S test shows no significant difference between the distributions. 
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lo 9 L 0 _5-8 kev (ergs s’) 

Fig. 5. The fraction of variable AGN (red circles) and galaxies (black squares) vs. X-ray 
luminosity. Error bars show the binomial errors at 90% significance level. The error bar on 
luminosity represents the bin size. The points are slightly offset for clarity. 
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Fig. 6. /f-barid magnitude vs. full-band X-ray flux for sources in the main CDF-S catalog 
of Xue et al. (2011), where filled symbols mark sources with significant variability (P x i < 
0.05). Sources classified as AGN, galaxies, and stars are marked with red, black, and blue 
symbols, respectively. Constant flux ratios are marked with diagonal lines: from top to 
bottom, \og(F x /Fn) — +1, — 1, and -2. The vertical dotted line marks the on-axis full-band 
flux limit for 20 net counts. The sources with ft-band magnitude lower limits at the top of 
the plot (s=»25% of CDF-S sources) have no /f-band counterpart. 
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Fig. 7. The distribution of max-to-min ratios for simulated light curves based on CDF-S 
galaxies (black histogram) and for measured light curves of variable galaxies (red histogram). 
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Fig. 8. The distribution of the bias-corrected normalized excess variance ojf X5 , X) , T for 
CDF-S sources classified as variable galaxies (red histogram), with non-variable galaxies 
(blue histogram) shown for reference. As expected, the distribution of non-variable galaxies, 
though affected by scatter due to statistical fluctuations, centers around zero (i.e., variability 
strength is attributable to Poisson fluctuations). The red dashed line marks the maximum 
variability expected from an XEB population. 
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Fig. 9. Postage-stamp images from the GOODS-S/GEMS HST V606 band for 20 variable 
galaxies. The label at the top of each image gives the source name. The numbers at the 
bottom of each image indicate the source number (X1D) in the main CDF-S catalog, the 
adopted redshift, and the logarithm of the full-band luminosity as calculated in §2. The 
circle overplotted on each image has a radius of 1.5 times the X-ray positional error, which is 
calculated at the 90% significance level, to illustrate whether an X-ray source is considered 
to be off-nuclear (Lehmer et al. 2006). Each image is 8" x 8" with the position of the X-ray 


Fig. 10. Postage-stamp images from the GOODS-S/GEMS HST V606 band for 72 non- 
variable galaxies. The labels and error circles are the same as for Figure 9. 
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Fig. 11. The bias-corrected normalized excess variance (crjj xs<x , tt ) vs. X-ray luminosity 
(Lo. 5 — 8tev) for CDF-S sources classified as AGN (red circles) and variable galaxies (black 
circles and stars) in the 4 Ms main catalog, where sources with £o. 5 _ aie v>10 41 ergs s -1 
arc fitted with a weighted least-squares regression (dashed line); the shaded area shows the 
dispersion around the fitted line. Black stars highlight six variable galaxies with significant 
XRB contribution to the total luminosity, discussed further in §6.2. Error bars include both 
measurement errors and errors associated with red-noise scatter and sampling (§4.1). Large 
blue squares mark the weighted means for each luminosity bin; error bars are the standard 
deviation of the data in each bin. Ultraluminous X-ray sources (ULX) are generally not 
expected above Lo. 5 - 8 tev ~ 10 41 ergs s - ' (vertical dotted line). 
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Fig. 12. The observed 2 10 keV X-ray power (Lj-iokev, tot) vs. that expected from an XRB 
population (La -lokcv.xmi) for variable (red stars) and non-variable (green squares) galaxies. 
Arrows designate upper limits on CDF-S AGN (open orange circles), 32 local 

galaxies (Colbert et al. 2004), and 20 local luminous infrared galaxies (Lehmer et al. 2010) 
are plotted for reference. Local galaxies lie on L2-i««v,tot L 2 . lokev.XRB (red solid line), 
with a measured dispersion of 0.34 (red dotted lines Lehmer et al. 2010). 
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Fig. 13. The bias-corrected normalized excess variance vs. accretion rate nor- 

malized by the Eddington rate (m= Ml Mem) for CDF-S sources classified as AGX (red 
circles) and variable galaxies (black circles). A weighted least-squares regression is fitted 
to all variable sources (dashed line), and the dotted lines show the dispersion around the 
best-fit line. 
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Fig. 14. The bias-corrected normalized excess variance (cr^ xscon ) vs. X-ray luminosity 
(Z-o. 5 — 8kcv) as in Figure 11 for sources with 2 < 1. Model relations are shown for a universal 
broken power-law PSD function, where the break frequency depends on Man and in. The 
relations are calculated for m ~ 9 x 1G~ 3 (dashed line) and rh * 4 x 10 1 (dash dot ted 
line), where Man — 10 s 10 l 0 A/© (low to high L o. 5 ~gk C v)- 






